ON "PASCAL", CODE GENERATION, AND THE CDC 6000 COMPUTER 



BY 



NIKLAUS WIRTH 



STAN-CS-72-257 
FEBRUARY 1972 



COMPUTER SCIENCE DEPARTMENT 

School of Humanities and Sciences 
STANFORD UNIVERSITY 



/fk 


^jy Ni p^ 


'vN 


"#: 


***->. "'"'< 




m 




I tO j; 


!^F 






*b, * \" 


<,. - 




%.Qa 


■< .-.--»<» . 


X. 


<?4NIZEOJ5? 





On "PASCAI", Code Generation, and the CDC 6000 Computer 

by Niklaus Wirth 



Abstract ; 

"PAoCAL" is a general purpose programming language with character- 
istics similar to ALGOL 60, but "with an enriched set of program- and 
data surac luring facilities. It has been implemented on the CDC 6000 
computei . Thic paper discusses selected topics of code generation, 
in parti- ular the selection of instruction sequences to represent 
simple operations on arithmetic, Boolean, and powerset operands. 
Methods to implement recursive procedures are briefly described, ana 
it is hinted that the more sophisticated solutions are not necessarily 
also the best. The CDC 6000 architecture appears as a frequent source 
of pitfalls and nuisances, and its main trouble spots are scrutinized 
and discussed. 



The preparation of this paper was made possible by support from the 
National Science Foundation, Grant number GJ-992, IBM Corporation, 
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On "PASCAL", Code Generation, and the CDC 6000 Computer 



This set of notes has a dual purpose. It is on the one hand 
directed to the user of the PASCAL compiler system who would like to 
gain some insight into the machine code which is generated for various 
basic operations. It is even recommended that he study these notes 
carefully, "because their understanding may prevent him from certain 
pitfalls which are inherent in the use of the CDC 6000 computer [1] . 

On the other hand the notes may be of interest to compiler writers 
in general, because they point out some problems and dilemmas and our 
choices of solutions. It becomes apparent that the choice of the code 
to be venerated is crucial for a good compiler system, and that i"t is 



- ^.l J. J. <jhi 



The true purpose of a higher-level language is that it allows a 
programmer to conceive his algorithms in terms of some convenient 
abstractions . For instance, he is given the opportunity to think in 
terms of familiar notions of numbers, of relations, and of repetitions, 
instead of having to express his program in terms of bitstrings, 
arithmetic instructions, and transfers of control. However, these 
abstractions are only truly useful, if he can assume that his implemen- 
tation observes all the properties which are commonly attributed to 
these abstractions, or else if it automatically issues a warning . As 
an example, when dealing with numbers in a uigu- j.evej_ _i_anguage, one 
sbould l^H^ i~ o assume all tb^ co r *iiron axioms of arithmetic to bold, ^f 
course this is not possible, since computers can only represent finite 
ranges of values. So one expects to receive a warning, if an operation 
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has trespassed the limits imposed by the implementation and an 
operation generates a result not in accord -with the rules governing 
the abstraction. So the system is expected to provide an error 
indication, e.g. if an overflow occurs in an addition, if a value is being 
assigned which lies outside the specified range of values of variables, 
or if an array index is used which lies outside the defined limits. 

Unfortunately, such potential warnings require the execution of 
additional instructions, which in general is costly. As far as range 
checking is concerned, they can be requested to be generated by the 
compiler for run-time execution by enabling so-called options . (The 
A-option generates assignment range checks, the X-option index checks.) 
They are relatively costly, but may speed up the finding of logical 






As far as irregularities of the arithmetic are concerned, one has 
become used to receive these warning signals automatically from the 
hardware, particularly because they are easily generated by the hardware, 
whereas a solution to detect overflow by software is usually beyond any 
reasonably economical feasibility . Unfortunately, the CDC computer 
fails to satisfy even the most modest expectations in this respect, and 
the effort to provide a system with security in the above sense was 
therefore a series of constant frustrations. Equally disappointing are 
some of the "features" of its floating-point arithmetic instructions. 

One can go only a relatively short distance in trying to correct 
mistakes of the hardware by means of software; otherwise a system 
becomes ridiculously inefficient and will not be used by conscientious 
programmers who are willing to take the peculiarities of a hardware 
into account and guarantee safety of their algorithms by analytical 



rather than experimental means. And this would have been against the 
intentions of PASCAL. So all that can reasonably be done is to 
elucidate the shortcomings and limitations of the hardware that are 
still transparent through the "software cover", and to make the programmer 
fully aware of them. And this is the purpose of this note. 

It concerns itself with the simple operations of integer and real 
arithmetic, with Boolean operations and with powersets. The reader is 
supposed to be familiar with the CDC COMPASS notation. The operands 
are usually assumed to have been brought into the XI and X2 registers. 
(If they were loaded into other registers, a corresponding renumbering 
is necessary which is, however, irrelevant to the operation itself) . 
Registers XI - X5 are used as a stack for intermediate results, whereas 
XO is used exclusively as local work register. 

Section 6 deals with the topic of implementing recursive procedures 
and the addressing of local variables. Although the general techniques 
are well-known, analysis of possible solutions and their experimental 
comparison yielded some noteworthy results. It is shown that attempts 
to make full use of available hardware features such as base registers 
may not necessarily lead to an optimal performance. Again, the 
instruction set of the CDC computer is hardly optimal to implement 
mechanisms for recursive procedures. Conspicuously absent is a sub- 
routine jump instruction which leaves the code invariant (reentrant) . 



2. Integer Arithmetic 

Data of type integer or of subranges thereof are represented by 
fixed-point binary numbers. Addition and subtraction are represented 
by the 

IXi Xj + Xk 
instructions. Other operations are implemented by short sequences of 
instructions, as outlined below. 

2.1 Mult iplic at ion 

Due to a recent change of the hardware, fixed-point multiplication 
can be performed by a single 

DX1 X1*X2 
instruction. It should, however, be noted that this instruction is 
essentially a floating-point instruction, and yields incorrect answers 
for fixed-point operands with |x| > 2 . This can be regarded as an 
overflow condition which is, alas, neither trapped nor indicated by the 
computer. A "safe code", checking against all imposed limits of operand,' 
and result, is quite elaborate and uneconomical by any standards, and was 
therefore not implemented. 

If one of the operands is a constant C being representable as 
either 

(2, h, 8, 16 ...) 

0> 5, 6, 9, 10, 12 ...) 
m > n 

(7, lh, 15 ...) 



1. 


c =2 n 


2. 


c =2 m +2 n 


3. 


c =2 m -2 n 



then the compiler generates the following code for the multiplication 
of XI by c : 



1: 


LX1 


n 


multiply by 


2 n 


2,5: 


LX1 
BXO 
LX1 
1X1 


n 

XI 
m-n 
Xl+XO 







Again, overflow conditions are simply ignored. Case 5 yields only 
correct results, if |xi*2 m | < 2 59 . 

2.2 Division ( div ) 

Integer division is represented by the instruction sequence 

EX1 XI 1 

EX2 X2 T pack 

MX2 X2 J 

FX1 X1/X2 divide 

UX1 B?,X1 

LX1 B7,X1 

m xi + xo J su PP ress ne §- zero 

and suffers from the same basic shortcoming as multiplication: an 

operand |x| > 2 yields an incorrect result. 

If the divisor is a constant c = 2 , the compiler again produces 

an "optimized" code, performing division by shifting. Unfortunately, 

a single right shift instruction is unsatisfactory, because it may 

generate a "negative" zero as result. Negative zeroes, however, must 

not be allowed to occur, since comparisons may yield wrong answers if 

applied to them. Thus, the optimized division is implemented as 

AX1 n divide by 2 n 

BXO XO-XO suppress 

1X1 Xl+XO negative zero 
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'ote that the unconditional generation and addition of a zero can be 



accomplished with a code that is not only shorter than a conditional 
jump, such as 

AVI v. 

*WJ- 11 

NZ XI, L 
SX1 BO 

I i ■ • * 

but also avoids the insertion of padding instructions (NOPs) for word 
boundary alignment. 

The D -opt ion provides an additional security measure against 
division by zero. It causes the compiler to insert a 

ZR X2, error 
jump instruction preceding every division instruction. (This applies 
to the Modulus operation as well.) It is particularly recommended in 
the case of integer division, where the actual divide instruction 
generates a "floating-point infinity" value, which is incorrectly 
treated by the subsequent conversion instructions and thereby represents 
a senseless result. 

2.J Modulus ( mod ) 

The modulus or remainder operation is defined as 
x mod y = x - (x div y) * y 
As it involves integer multiplication and division operations, it 

suffers again from the same deficiencies of the 6000 arithmetic. Its 

corresponding code is: 

PXD XI 

PX6 X2 

nx6 x6 

FX6 X0/X6 

UXt- B7,X6 

1x6 B7,x6 

DXb XD*X2 

1x1 xi-x6 



2.h Sign inversion 

The use of one's complement representation for negative numbers 

BX1 -XI 

unsatisfactory, because it might generate a "negative zero". So we use 

BXO XO-XO 
1X1 X0-X1 

2.5 Comparisons 

Since the computer does not offer a compare instruction, subtraction 
has to be used; this has primarily the disadvantage of generating wrong 
results in the case of overflow. The cases of testing for equality and 
inequality are handled correctly, because the one's complement addition 
generates an end-arcund-carry in the case of "negative overflow", thus 
maintaining a result indicating inequality. Note that the Boolean 
subtraction 

BX1 X1-X2 
cannot be used, because a comparison of xl and x2 = -xl would yield 
a zero result, thus indicating equality. 

Whereas equality testing is "safe" with the 

1X1 X1-X2 
instruction ignoring overflows, this is not the case for the tests of 

ordering (xl < x2) by subtraction and subsequent inspection of the 

1 1 59 

sign bit. The reason is that if overflow occurs, i.e., |xl-x2| > 2." , 

then the sign bit will be the opposite of the true sign. This situation 

is quite hopeless, since overflow is in no simple way detectable on this 

machine. In order to obtain a (sign) bit representing the relation 



x < y for any values x,y , the following algorithm can be used: 

1. Compare the signs of the two operands. 

2s If the - "' are different- then the result is obvious. 

3. If they are equal, the subtraction x-y can be performed 

without danger of overflow, and x-y < is the result . 

A minimal instruction sequence to perform these operations and avoiding 

the use of undesirable jump instructions is 

BXO X1-X2 compare sign bits 

1X2 X1-X2 

BX1 X0*X1 if unequal, choose sign of XI 

BX2 -X0*X2 if equal, choose sign of X1-X2 

BX1 X1+X2 

Now the sign bit of XI is 1 , if XI < X2 , and otherwise. Still, 

the effort to perform a faultless comparison is formidably cumbersome, 

and the PASCAL compiler does not generate it. The programmer is left 

with the responsibility to verify that for every comparison of x and 

|x-y| < 2 59 . 

2.6 Taking the absolute value ( ABS ) 

The code used to take an absolute value is designed to avoid jump 

instructions, not only because they are long and slow, but because they 

usually introduce NOP instructions for alignment. 

BXO XI 

AXC 59 generate 60 sign bits 

BX1 X0-X1 

2.7 Testing for even or odd (ODD) 

Since one's complement representation is used for negative numbers, 
the least significant bit of the operand must be compared with its sign 
bit: 
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BXO XI 
LXO 59 
BX1 XI-XO 

This leaves the sign-bit of XI equal to 1 , if XI was odd, and 

otherwise. 

The compiler "optimizes" in the case of ODD(x) with x being of 
a subrange type with only non-negative values. It then generates the 
single instruction 
LX1 59 

2.8 Summary 

The foregoing explanations reveal that the absence of any overflow 
indication makes analytical verifications necessary that guarantee the 
non- occurrence of these conditions. An effective aid in experimental 
testing is the A-option, causing interval check instructions to be 
generated with every assignment to a variable that is declared to be of 
a subrange type. The A-option is activated by the "comment" 

[$A+ . . . } 
and causes the code for an assignment to a variable 

VAR V: a..b 
to become: 

SX7 * location identification for error trap 

SX1 a 

1X0 X6-X1 

SX1 b 

1X1 X1-X6 

EXO X1+X0 

NG XO, error jump to error routine 

SA6 V 

It should be noticed that unfortunately the attractive and shorter code 

sequence 



SX7 * 

SXO x6-a 

SX1 X6-b-l 

BXO -Xl+XO 

NG XO, error 

SA6 V ' 

cannot be used., because the instructions 

SXi X j + K 

perform an 18-bit arithmetic ignoring the leading 1+2 bits of the 

register Xj which --of course -- is not in the spirit of a check. 

This ignoring rather than checking of the leading bits in 18-bit 

arithmetic is the reason why the so-called "increment" instructions 

cannot be used by the PASCAL compiler, except in the following special 

circumstance: if a variable x is declared of a subrange whose limits 

17 
are both less than 2 in absolute value, then the assignment statement 

x := x+ k 

is compiled as 

SA1 x 
SX6 Xl+k 
SA6 x 
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3. Floating-point Arithmetic 

The PASCAL compiler uses the complete set of F-instructions for 
arithmetic with values of type "real". Comparison is performed by- 
subtraction due to the lack of a compare instruction. This is possible 
without handicap since the occurrence of overflow generates a signed 
"infinity" -value, but no immediate trap. Sign inversion is represented 

by 

BXO XO-XO generate zero 
DO. X0-X1 

and the absolute value function by 

BXO XI 

axo 59 

BXl X0-X1 
Arithmetic with the F-instruction possesses some peculiar properties 

which will briefly" be reviewed, and has for instance the consequence that 

-.-_t- "i -i-,.-.r V1.-.+- ^- ~ -> .-. r. ,- .-. >- • •> ■, • .„.,- -> ,.- .. , -• r '■ i -j ^-•■ n f t vir>"- • 

computed by an Finstruction. The trouble arises from the fact 

that F-arithmetic truncates without rounding, and F-addition truncates 

without post -normalization. Every addition is therefore compiled into 

two instructions: 

FX1 X1+X2 add/ subtract 
KX1 XI post -normalize 

If the two values 

a - 1720 1*0... 00B - 1.0 

rt -U8 
b = 1717 17...77B = 1.0-2 

are compared by subtraction 

FXO X1-X2 a-b 

the result is 
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1720 to 00 / 00 

- 1720 37 77 / to 

1720 00 00 / 1+0 

"where the slash marks the separation between the lower and the upper 
half of the 96-bit accumulat or . The result is although the two 
operands were different. 

Notice that subtracting 0.5 from both a and b , and then 
computing their difference, yields 

a - 0-5 = 0-5 : 1717 ho 00 

_ii_8 

b - 0.5 = 0.5-2 : 1716 77 76 

1717 to 00/00 

- 1717 37 77/00 

171 7 OP.*. ••.... ••■•P 1 / 00 ° 

i.e. j a difference which is not zero. Thus the result does not only 
depend on the true result, but also on the values of the operands. 
This unpleasant property of the CDC F-arithmetic stems from the fact 
that automatic post -normalization is absent. 

3.1 Rounding 

It was at one time hoped that this defect could be avoided by 
letting the PASCAL compiler automatically generate R- instruct ions, 
which include a certain kind of rounding. However, R-arithmetic 
turned out to feature some even stranger properties, so that it was 
decided not to use R- in struct ions. In order to point these features 
cut, a brief review over R-arithmetic is necessary: 
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The R-instructions differ from the F- instruct ions only insofar 
as a 1-bit is appended to normalized operands before the arithmetic 

operation is performed. Thus for instance the subtraction of 

_ii.8 

b = 1.0-2 from a = 1.0 yields 

4 
1720 1*0 00 / ko 

-1720 37 U/60 




1-bits appended 



1720 00 00 / 60 

which of course is still zero. 

The principal defect with "CDC -rounding", however, is that its 
effect is unpredictably either the addition of l/2 or l/k in the 
last position, because rounding takes place before instead of after 
normalization (which must again be performed by a separate instruction) 

five-bit number representation: 

16 = 10000 / 1«— _ inserted 

+17 = 10001 / 1 < — round -bits 

33 = 100010 / 

10001 / = 3^ 



31 = mil / 1 4 inserted round -bit 

+2 = 00010 / 

33 = 100001 / 1 

10000 / 1 = 32 

In the first case, the pre-rounding results in correct rounding of the 

■pre-roundinp has no effect . 

The same phenomenon can be observed in the cases of multiplication 
and division. The following example again uses a five-bit number 

representat ion : 
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round -bit 



15x12 = 11110 / 1 x 11000 



11110 / 1«— 

+ 01111 / Qlf. 



101101 / 11 
10110 / 111 



176 



18 X 10 



round-bit 



10010 / 1 x 


10100 


10010 / 1« 

1 00 / 101 £ 


1 











10111 / 001 



l8iv 



In the first case, the roundin^ effect is nil, leaving the inexactly 
representable value ISO be an unrounded Y[(j\ in the latter case the 
rounding effect transforms 180 into the value 18^. (Suitable adjust- 
ment of exponents is not shown here.) 

A method introducing proper rounding instead of "CDC -rounding" 
relies on the use of the D-instruction set [2]. Whereas the F-instructions 
yield the high-order kQ bits of the 96-bit accumulator, the D- instruction;; 
yield the low-order hQ bits with a suitably adjusted exponent, thereby 
allowing access to a double precision result. 

Notice that it is an ingeniously efficient method to compute a 
double precision result by 

"1 f-i, .-.wi-t-v'4- tv>t 4-Vk rt T\"D ■v»^> mi "i -J- qm^ A -1 c"-r-\/-\ o £± n-P 4-Vi/^ ~\ <*\t»t Vi q T ~P 

J_ • ^\->IlLjJU.\J±.ll£± i OllC ±JJl —± COUJ-U Ctliu. u.a.u_ywuC wj. UULi _i_v>W .law,.!— i- 

( "p_ instruct"? on ) - fb^p 

2. computing the same again and dispose of the high half 
(D-instruction) - 
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This computer allows it to be done in no other way! 

The PASCAL compiler will generate the following code for floating- 
point operations, depending on the choice of the R- opt ion : 

ON 



R-option 




OFF 


x + y 


FX1 
Ml 


XI + X2 
XI 


x * y 


FX1 


XI * X2 


x / y 


m 


XI / X2 





FXO 


XI + X2 




NXO 


XO 




DX1 


XI + X2 




RX1 


XI + XO 




MX1 


XI 




FXO 


XI * X2 




DX1 


XI * X2 




RX1 


XI + XO 




RX1 


XI / X2 



Examples of addition / subtraction: 



1. 



1 '° -48 
1.0-2 


1720 
1717 
1720 


40. .. 
77... 

37.-. 


. . .00 / 00. 

...77 / 00. 

...77 / 40- 


* • • • • \J 



• • • • • \J 


= 


1720 
1640 


00... 
40... 


...00 / 40. 
...00 


• • • » • \J 

after addition of high 
and low 



-48 
2. Take a = 1.0 and b = 2 , then subtract a-b : 



F -subtract ion yields 

a = 1720 40 00/00... 

b = 1720 00 00 / UP... 








1720 37 ■ 



...77 / 40... 



which, after normalization, is 
1717 77 76 



1.0-2 



■hi 



R-subtraction inserts a 1-bit after the slash in the first operand, 
ana thus yields "the result 

1720 40 00 = 1.0 exactly 
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The combined use of F and D instructions yields the true result, 
because the normalization instruction left shifts the high order 
result to 

1717 77 76 

whereafter a "rounded" - addition is used to add the correction 
+ 1717 00 01 

yielding 

-J+8 
1717 77 77 = 1.0-2 

3.2 Conversion from fixed to floating-point (integer to real) 

Wherever a real operand is permissible, PASCAL allows the specifi- 
cation of an operand of type integer as well. However, the compiler is 
then forced to generate the necessary representation conversion instruc- 
tions, which are not only time-consuming, but potentially hazardous. 
It is therefore recommended to avoid "mixed-mode" arithmetic expressions 
wherever possible. The generated conversion instructions are 

PX1 BO, XI pack with zero exponent 
UX1 BO, XI normalize 

The result of this conversion is wrong, whenever the integer operand 

in XI is larger or equal to 2 in absolute value, since the 

exponent bits are simply ignored by the P instruction. A test to verify 

that the operand is within bounds could be compiled as 



BXO 


XI 


AXO 


k8 


NZ 


X0, error 


3,S J-.I.V 


-"i oyi -f- y~i *V> r\ 1 



>vi C\ r* .'> c 



■+- *1 ir 4-Vi f^v» 4- Vl O AAriirpyici -i nri "i T O Q 1 "P 
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3-3 Conversion from floating to fixed-point (real to integer) 

PASCAL does not provide for any implicit real to integer conversion. 

However, the standard function TKUNC(x) allows to truncate the 

fractional part of a real number. The used code is: 

UX1 B7,X1 

LX1 B7,X1 

BXO XO-XO "V avoid 

1X1 X1+X0 J negative zero 

The result of this conversion is again wrong, if |x | > 2 . 
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h. Boolean Operations 

The standard type Boolean is defined in PASCAL as 
type Boolean = (false, true) 
Since the values of all scalar types are mapped onto the integers 
0,1, 2, . . . , the values false and true are represented by the numbers 
and 1 respectively. 

The operations A and v are implemented by the Boolean AND and 

OR instruction, namely 

BX1 X1*X2 and 
BX1 X1+X2 

Negat ion is performed by 

MXO 59 
BX1 -X0-X1 

If a relation has to be assigned to a Boolean variable, e.g. 

b : = x < y 

then a sequence of instructions is necessary to obtain a or 1 

value. Again every effort is made to avoid the use of jumps. The 

following code is used in the above assignment; leaving a Boolean value 

in XI . 

FX1 X1-X2 x-y 

MXO 1 

BX1 X0*X1 Extract sign bit 

LX1 1 move it to correct position 

Analogous code is generated for the relations > , < , and > . But 

unfortunately the equality relations cannot be reasonably implemented 

without a jump; in the assignment 

b := x = y 
the following instructions are generated: 
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FXO 


X1-X2 


BX1 


XO-XO 


NZ 


XI, L 


SX1 


1 


J_J • • • 





Boolean comparisons, although occurring rather infrequently, are treated 

as special cases, because a simpler and shorter code is applicable: 

p < q BX1 -X1*X2 

p < q BX1 -X2*X1 

MXO 59 I j.. 
BX1 -X0-X1 )™^™ 

p £ q BX1 X1-X2 

The remaining three relations are compiled analogously. 
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5* Powerset Operations 

PASCAL 6000 restricts powerset types to be built only on base sets 
with less than 59 components. This allows a powerset value S to be 
represented by one "word", in which the i-th bit indicates the presence 
(1) or absence (0) of the element i in S . 

5-1 Generation of the Singleton Set [i] 

Assume that i is loaded into register XI, then 

...SSI. XI. 
SX1 1 
LX1 B7,X1 

Notice that the numbering of bits starts with at the low order end. 

This choice was made in order to be able to load powerset constants 

with small valued components (less than 18) by a single SXi instruction. 

5«2 Set Intersection, Union, and Difference 

These three operations are implemented by a single instruction 

intersection BX1 X1*X2 
union BX1 X1+X2 
difference BX1 -X2*X1 

5-5 Set Membership (in) 

The relation i in S is implemented by shifting the bit representing 
i into the sign position which can be tested: 



SB7 


XI 


i 


* in 

fiA-L 


B7,X2 


S 


LX1 


59 





11 the expression i is in the form of a constant c , then the compiler 
generates of course only the single instruction 
LX1 59-c 
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5.k Set Comparison 

Sets can be compared for equality and inclusion. Equality is 
tested by a Boolean subtraction 

BXO X1-X2 
and a subsequent zero test. Note that the peculiar property of the 
zero test to recognize a word with either 60 zero-bits or 60 one-bits 
as a zero is responsible for the restriction that powersets may contain 
at most 59 instead of 60 elements. If sets with 60 components were 
allowed, then a full set and an empty set would not be distinguishable 
by a single subtraction followed by a zero-test. 

Inclusion expressed as x < y and meaning x c y , is implemented 
by the single instruction 

BXO -X1*X2 
which is followed by a zero-test instruction. The same instruction is 
used for the relation x > y , whereas strict inclusion (x c y) is not 
implemented. 
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Some Exercises Addressed to the CDC 6000 Expert 

1. Is the following code to represent the function trunc(Xl) 

acceptable? If not, why? 

BXO XO-XO 

FXO XO 

FX1 X1+X0 

UX1 B7,X1 

NZ B7, overflow 

2. Is the following code for XI mod X2 acceptable? If so, prove it, 

FX1 XI 

FX2 X2 

NX6 X2 

FX6 X1/X6 

BXO XO 

FX6 X6+X0 

DX6 x6*X2 

fx6 xi-x6 

uxi x6 



3. Why can the instructions 

BXO X1-X2 

ZR XO, equal 

not be used to represent a comparison XI = X2 ? Prove that 

1X0 X1-X2 

ZR XO, equal 



always yields the correct action. 
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(•>. Implementation of Recursive Procedures 

The language PASCAL has been carefully designed so that dynamic 
storage allocation is not required, with the following two exceptions: 

1. Variables local to procedures may be allocated storage only 
when the procedure is called, and 

2. Components of class variables are allocated storage by calling 
the standard procedure "alloc". An area of store is allocated 
to the entire class variable as soon as the procedure is 
called to which the class is local. 

In this section we will briefly review the well-known techniques for 
handling recursive procedure calls and of allocating storage to their 
local quantities, and discuss the code selected to represent the 
procedure call mechanism. 

Due to the first -in last-out nature of the hierarchy of activated 
procedures a stack may be used to allocate local variables. This is of 
great advantage, since storage retrieval is trivial in the case of 
stacks, resulting in low storage management overhead. We consider the 
set of local variables of each activated procedure as a record (often 
called "data segment") in the stack. Since their lengths may all be 
different, the most convenient method to thread the way back through 
such a stack is by constructing a chain of pointers linking the records. 
Every record then contains a "header" containing 

1. the link to the previous record, and 

2. the (frozen) program status (counter) of the calling procedure. 

variables are addressed relative to the origin of the record of 
which they are a part. The origin address is unknown at compile-time, 



and must be determined at run-time. This can be done by descending 
through the link chain, until the desired record is reached. But how 
is the desired record recognized? The most straight -forward method 
which interprets the scope rules of an ALGOL block structure correctly 
is probably the following: 

Method I : 

1. Define the level of an object to be 1 greater than the 
level of the procedure to which it is local. The level of 
the main program is . 

2. Indicate the level of each record (equal to the level of its 
components) in its heading. 

5. Whenever an object on level i has to be accessed, the record 
containing it is found by descending down the chain of links 
until the first occurrence of a level indicator with value i 
is found. 
This accessing method has the obvious drawback of inefficiency (and of 
not being applicable in the case of parametric procedures) . A slight 
modification, however, improves efficiency and generalizes to parametric 
procedures. 

Method II ; 

Instead of indicating levels explicitly in the record headings, 
a second link chain is constructed connnecting each record A with its 
static ancestor, i.e., with the record B of the procedure in which A 
was declared locally. In order to distinguish the two link chains, the 
former is called the "dynamic link" and the latter the "static link". 
An example of a state of computation is shown below for a given -- 

admittedly not very realistic -- program. 
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var vO; 

procedure QO ( procedure X) ; 
var wl; 
procedure Ql; 

var w2; 
begin w2 := wl+vO; X 
end ; 
begin wl := vO; Ql 
end ; 

procedure PO; 
var vl; 
procedure PI; 

var v2; 
begin v2 := vl+vO; QO(PO) 
end ; 
begin vl := vO; PI 
end ; 
begin {main program} vO := 0; PO 
end. 




stack 



Ql 



QO 



PI 



PO 



Ql 



QO 



PI 



PO 



Ql 



Qw 



PI 



PO 



main 



dynamic 
link 



:) 



'S 



■ 



') 



i 



• 



= N 



} 



static 

link 



Method III : 

Although the use of a high-speed index register to represent the 
origin of the link chains improves access speed significantly, the 
process of descending down the static chain to the record (data segment) 
with the desired level is relatively time-consuming. An ingenious 
device to reduce access time was introduced by Dijkstra [k] and is now 
widely used in compilers for block-structured languages. The device 



2i> 



is an array of base addresses, called the Display D , which is at any 
time a copy of the static chain. If an object at level i is to be 
accessed, the origin address of its data segment is quickly obtained 
as D. . The method is particularly attractive for computers with a 
set of high-speed index registers which can be used as the Display. 
The price for this increase in access speed -- apart from the reservation 
of registers -- is the setting and updating of the Display each time a 
procedure is called and terminated. To be more specific, the necessary 
actions are as follows: 

1. if an actual procedure of level i is called, D. has to be 
set; 

2 . if control is returned from a procedure at level i to one 
at level j , (j > i) , D. ...D. have to be reassigned; 

J. if a formal procedure at level i is called from a procedure 
at level j , D. ...D, have to be reassigned, where k is 
the level on which the static link emerging from the calling 
and the called procedures merge. Since k is not known at 
the time the procedure declaration is compiled, k can be 
chosen as zero without significant loss in efficiency. 
This scheme was used in the implementation of PASCAL 6000. It is 
described in Reference h. Registers Bl . . . B5 are used as the Display, 
B5 is the origin of the link chains, and B6 is the pointer to the top 
of the stack. The compiled instructions are the following: 
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Procedure call of P 



SX7 



L 
P 



L SBj 



B5 



SAl B5 "1 
SB(j-l) XI J 



SAl 
SBi 



XI 
XI 



-\ 



J 



save return address 
and jump 



\ update the display, if j > i 



Procedure entry: 



SBi 


XI 


SAl 
SB(i 


XI 1 
-i) XI J 


SAl 
SBI 


XI I 
XI J 


SA7 
SX7 
SAT 


B5+-2 

B(i-l) 

B6 


SX7 
SA7 
SBi 

SB5 

SB6 


B5 

B6+1 

b6 

B6 

B6+L 



} 



prolog, entry for calls 
of formal procedures 



update display 

save return address in header 
save static link 

save dynamic link 

new display entry 

T 

top of stack, L = data segment length 



Procedure exit: 



SB6 


B5 


SAl 


B5+1 


SB5 


XI 


SAl 


B5+2 


SB7 


XI 


JP 


B7+0 



} 



reset top of stack 
reset T 

fetch return address and jump 



Notice that global variables in the main program are assigned absolute 
addresses. Since BO = , they can be considered as based on BO . 
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In the first half of 1971, Prof. C. A. R. Hoare and his collaborators 
modified and bootstrapped the PASCAL compiler for the ICL 1966 computer [6], 
une u± one more ssi^iij-xxuaiiL. aioeiaoxuHo <^uii>-cj.iicvj. uuc gjajuiuuuiuu ^j. 
the Display, due to the fact that the ICL computer has no set of index 
registers that are available for a Display, and since the use of a 
Display was not considered to be an advantage, in this case. During 
a visit of Prof. Hoare in July 1971, he suggested that maybe even with 
a register set available for the Display, the benefits gained should be 
investigated. His suggestion was certainly valid, since variables 
either global or local to the most recently called procedure could be 
accessed with the same speed even without a Display. Thus the gain from 
a Display is limited to faster access of objects at intermediate levels, 
while the price is the updating at every call regardless of whether such 
objects are accessed or not. A superficial look at the PASCAL compiler 
itself showed that accesses to such intermediate level objects were 
indeed relatively rare, and it was decided to generate a version that 
would not use a Display (Method II) . This version still uses the address 
register B5 as origin of the link chains (and base address of the most 
local data segment) and B6 as pointer to the top of the stack. The 
generated code is: 

Procedure call P 

* X6 := base of environment of P 
SX7 L 

m P 

L 
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Procedure entry: 

SXO B5 

LXO 18 



-\ 



pack and store 
BX7 X7+X0 ) dynamic link and 



SB5 B6 

SA7 B5+1 -> 

SB6 B7+L 

* SA6 B5 



return address 

stack pointer 
static link 



Procedure exit: 



SA1 


B5+1 


sb6 


B5 


SB7 


XI 


1X1 


he 


SB 5 


XI 


JP 


B7+0 



~\ 



fetch and unpack 
dynamic link and 
return address 



Fetching an object x at level j from code at level i : 

1) J=0: SA1 BO+x 

2) j = i: SA1 B5+x 



3) < j < i: 



SA1 B5 
SA1 XI 
SA1 Xl+x 



repeated i-j-1 times 



A comparison of the codes generated by the two compilers shows that gains 
and losses of execution speed should be measured, but also those of code 
length. The shorter codes for procedure entry (2 - 2-j- words vs. 
h - 6 words), procedure exit (2 vs. 3 words), and procedure calls 
(no updating of display) are very attractive, particularly in a compiler 
where space is more on a premium than time, (it should be noted that 
the instructions marked with an asterisk can be omitted in the call or 
the entry code of procedures declared on the first level) . Of course 
it must be kept in mind that the decision about which compiler is to be 
preferred depends not only on the weighting of space vs. time, but even 
more on the programs to be processed. But it is obvious that if the 
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majority of these programs rarely use nested, procedure declarations, 
and often call procedures on the same level, then the compiler without 
Display is to be preferred. The compiler itself, although featuring 
nested procedure declarations, but seldom accessing intermediate level 
variables, belongs to this class. Comparisons of code generated by the 
two compiler versions produced the following results: 

1. The efficiency of codes not using a Display is in the average 
slightly higher (the compiler itself runs about 1.5$ faster). 

2. The size of codes not using a Display is smaller (by about 
hi, measured on 25 sample programs, about 6<f in the case of 
the compiler's code). 

3. The compiler program itself is slightly less complex without 
Display. 

This episode where a more sophisticated method was abandoned in favor 
of a simpler and more direct technique could well be added to the list 
of D. Khuth's examples of adverse influences of "computer science" on 
"computer usage" [5]« Their common characteristic is that improved 
methods are adopted without closer inspection of the nature and direction 
of the improvement, and without analysis of the circumstances to be 
improved. An interesting fact is that the Burroughs B5500 computer — 
specifically designed for ALGOL implementation -- did contain exactly 
the two base registers required to efficiently address objects at 
levels and i . Unfortunately, addressing of intermediate level 
objects was impossible due to the software; this deficiency was 
justifiably criticized. The remedy adopted in the successor BopOG was, 
however, not a correction of the deficient software, but the inclusion 
of a full set of high-speed registers to serve as Display. 
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7. Summary of the Main Trouble Spots of the CDC 6000 Architecture 

1. Use of one's complement arithmetic. In order to keep comparisons 
simple and efficient, the occurrence of negative zeroes must be 
prohibited. (Note that PL and NG test the sign bit only.) Various 
optimizations are more cumbersome and less effective, because 
negative zeroes must be suppressed by additional instructions. 
Some instructions are themselves unsafe against the generation 

of -0 ! 

2. No overflow check on fixed-point arithmetic. This lack is very 
serious and may cause wrong restuls in totally unexpected 
situations. Overflow check by software is prohibitive. 

3. No compare instructions. The use of subtraction may cause wrong 
results, unless expensive precautions are taken. 

h. Use of ^8-bit multiplier and divider for fixed-point 60-bit numbers 
without warning of possible "overflow" of operands. 

5. Floating-point addition and subtraction without automatic post- 
normalization . 

6. Floating-point arithmetic with rounding of operands instead of 
rounding postnormalized results. 

7. No subroutine jump instruction depositing the program counter P 

in an operand register, and no return jump loading P from a general 
operand register. This defect requires the use of 3 instructions 
each to jump and deposit a return address, and to retrieve it and 
return, whereas many other computers need only a single instruction 
for these purposes. 



Conc lusions 

When considering these complaints., the reader should bear in mind 
that this computer's architecture was conceived in the very early 1960*s. 
The CDC 6600 machine was a very advanced design for a special purpose: 
fast number crunching. The design relied heavily on the use of several 
arithmetic units working simultaneously ("in parallel") . Integer 
arithmetic was considered as almost dispensible, and overflow interrupts 
as undesirable, because of the impossibility to mirror the present state 
of the entire machine by a simple program counter and of resuming compu- 
tation. The use of simultaneously operating units is apparently also 
made responsible for the otherwise incomprehensible absence of post- 
normalization, namely because the unit for floaint-point addition does 
no u con^axn a leio—Si-iixTi cxrcuxury. n j.ew years a-ater, one v^jjG w^r^o 
(and 65OO) computers were announced; they were to have the same instruction 
set as the 6600, but only one conventional integrated arithmetic -logical 
unit. Although the "reasons" for the absence of interrupts and post- 
normalization had vanished, these "features" were retained in the n ame 
of compatibility . It was apparently considered most important that 
pitfall loaded programs could be transported to the new machines at no 
extra cost. This policy of staying "upward compatible with all previous 
mistakes" was sternly maintained when the successor to the 6000 series 
was announced in 1971' 

This attitude, which is by no means atypical among computer 
manufacturers, makes it doubtful whether any progress toward more 
reliable and more efficient computing will ever be achievable. It 
does not seem so, until the computer consumers' attitudes will no longer 

-ill 0+ -1 -FSr -f- Vi o -r»-r»ac" a-trf- mown -Pa ri+iTvovQ ' t*/~*T -i r> -1 a a TVi enr -i n "hnvm t*t*i ~\ 1 nrvf* 
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change before they are made aware of the hidden cost involved in using 
the present equipment . I am convinced that the cost incurred by the 
programmers having to discover bugs the hard way by reprogramming 
repeatedly, and having to reexecute programs many times until they 
were believed to be correct, is incomparably higher than the reduction 
in cost due to staying compatible with outdated architectures. The 
project to develop the PASCAL compiler for the CDC 6000 computer 
unfortunately provided ample support for this conviction. 
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005001 
005001 
00500** 
005007 
005010 
005012 
005076 
005105 
005112 
005117 
005117 
005127 
005127 
005140 
005145 
fr05t54 
005154 
005160 
005171 
005200 
005210 
005212 
005212 
005216 
05223 
05223 

005074 



005075 
005076 

005077 
005100 

005101 
005102 



t$C + Til EXPRESSIONS AND ASSIGNMENTS > 
VAR I,J,K1 INTEGER} 
X,Y,Z« REAL? 
Nl 0..9999J 

P,Qi boolean; 

BEGIN i REAL ARITHMETIC > 

X 1= 1.0J Y 1= X + 3.14159; Z 

x i= x + cy + <z + (i.o + x) >>; 

X 1= ABS{+Y)J Y 1= SQR<X>; Z 1= 
C$R+ ROUNOEO REAL ARITHMETIC} 

x i= y + z; x i= y*z; x i= y/z; 

i INTEGER ARITHMETIC > 

t= i; j i= I + 100; k i= i * j; 
t= (-j) mod k; j i= SQR(J); 
i= trunc(X) ; z i= ij x i= i/j; 
boolean arithmetic > 



Is X»Y ♦ X/Y 



x; 



I 

K 

I 
t 
p 
p 
p 
p 

Q 

i 

I 
J 

N 
END 



SA7 
SX7 
SA7 



K 1= I DIV J! 



<Qvp>; 



BX6 
SA6 
NO 



SA1 
SA2 

FX1 

NX6 
SA6 



SA1 
SA2 

FX1 
SA2 
NO 



j; 
j; 
j; 



Q 
Q 
Q 



i= true; q. i= p a ■ 
is x = y; p t= i = 
i= x < y; p i= i < 

1= X < YJ P 1= I < 

1= ODD(I> ; 

OPTIMIZATION OF INTEGER 

1= 1*8 ♦ j»io; 

1= I OIV 8 - N OIV 2? K 
1= I ♦ 100 



85*80 
B5+80 
80+005000 



1= P s 
Is P < 

1= p < 



q; 

q; 
q; 



ARITHMETIC > 



1= I MOD 16; 



SB6 B5+000001 
SA1 80+005225 



XI 
80+005004 



80+005004 
30+005226 

X1+X2 
80, XI 
80+005005 



005103 

005104 

005105 
005106 
005107 

005110 



SA3 
FX2 
FX1 

NX6 
SA6 
NO 



80+005004 
BO+00 5005 

X1*X2 
BO+005004 



005111 



80+005005 

X2/X3 

X1+X2 

BO, XI 
BO+005006 



SA1 BO+005004 

SA2 80+005005 

SA3 80+005006 

SA4 BO+005225 

SA5 BO+005004 

FX4 X4+X5 

NX4 80, X4 

FX3 X3+X4 

NX3 B0,X3 

FX2 X2+X3 

NX2 80, X2 

FX1 X1+X2 

NX6 BO, XI 

SA6 BO+005004 
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005112 








SA1 


80 + 005005 




BXO 


XI 




AXQ 


fi 


005113 








BX6 


xo-xi 




SA6 


BO+00 50 4 


005114 


NO 










SA1 


80+005004 




FX6 


X1*X1 




NO 




005115 








SA6 


BO+005005 




SA1 


80+005004 


005116 








BXO 


XO-XO 




1X6 


XO-XI 


005117 


SA6 


B 0+005006 








SA1 


BO+005005 




SA2 


BO+005006 


005120 








FXfl 


X1+X2 




NXO 


B0,X0 




DX1 


X1+X2 




RX1 


X0+X1 


005121 








NX6 


BO, XI 




SA6 


BO+0050Q4 


005122 


NO 










SA1 


BO+005005 




SA2 


BO+005006 


005123 








FXO 


X1*X2 




DX1 


X1»X2 




RX6 


X0+X1 




NO 




00512<t 








SA6 


60+005004 




SA1 


BO+0050 05 


005125 








SA2 


BO+005006 




RX6 


X1/X2 




NO 




005126 








SA6 


80+005004 




SX6 


BO+000001 


005127 








SA6 


80+005001 





SA1 


BG+005001 


005130 








SXO 


B0+0001<t4 




1X6 


Xl+XO 




NO 




005131 








SA6 


B 0+005002 




SA1 


BO+005001 


005132 








SA2 


80+005002 




DX1 


X1+X2 




BXO 


XO-XO 


005133 








1X6 


Xl+XO 




SA6 


B0+Q05003 


005134 


NO 










SA1 


BO+005001 




SA2 


80 + 005002 


005135 








PX2 


B0,X2 




NX2 


80, X2 




PX1 


80, XI 




FX1 


X1/X2 


005136 








UX1 


87, XI 




LX1 


B7.X1 




8X0 


XO-XO 




1X6 


Xl+XO 


005137 








SA6 


80+005003 




SA1 


80+005002 


005140 








BXO 


XO-XO 




1X1 


XO-XI 




SA2 


80+005003 


005141 








PX6 


80, X2 




NX6 


B0,X6 




PXO 


BO, XI 




FX6 


xo/xe 


005142 








UX6 


87, X6 




LX6 


B7,X6 




0X6 


X2*X6 




1X6 


X1-X6 


005143 








SA6 


BO+005003 




SA1 


80+005002 


005i<*<» 








DX6 


X1*X1 




SA6 


B0+0C5002 




NO 
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005145 

005146 

005147 

005150 

005151 
005152 



005153 

005154 
005155 
005156 

005157 

005160 
005161 

005162 



SA1 
UX1 
LX1 

BXO 
1X6 

SA6 



SA1 
BX6 
PX6 

NX6 
SA6 
MO 



SA1 
SA2 

PX2 
NX2 
PX1 
NX1 

RX6 
SA6 
NO 



SX6 
SA6 



SA1 
SA2 

SA3 
8X2 
HXO 

8X2 

BX6 
SA6 



SA1 
SA2 

1X0 

MX6 
NZ 

SX6 



80+005004 
87, XI 
B7,X1 

XO-XO 
Xi+XO 
BO+0Q50Q1 



80+005001 

XI 

80, X6 

B0,X6 
80+005006 



B0+Q050Q1 
80+005002 

B0,X2 
80, X2 
BO, XI 
BO, XI 

X1/X2 
80+005004 



80+000001 
80+005010 



80+005010 
80+005011 

80+005010 

X2vX3 

73 

"•X2-X0 

X1aX2 

BO+005011 



80+005004 
BO+005005 

X1-X2 

00 

XO, 005163 

80+000001 





NO 




005163 








SA6 


80+005010 




SA1 


80+005001 


005164 








SA2 


BO+005002 




IXO 


X1-X2 




HX6 


00 


005165 








NZ 


XO, 005166 




SX6 


80+000001 


005166 








SA6 


80+005010 




SA1 


80+005010 


005167 








SA2 


80+005011 




BX1 


X1-X2 




HXO 


73 


005170 








8X6 


"Xl-XO 




SA6 


BO+005011 




NO 




005171 








SA1 


80+005004 




SA2 


80+005005 


005172 








FX1 


X1-X2 




HXO 


01 




BX6 


XOaXI 




LX6 


01 


005173 








SA6 


80+005010 




SA1 


80+005001 


005174 








SA2 


80+005002 




1X1 


X1-X2 




HXO 


01 


005175 








BX6 


X0*X1 




LX6 


01 


A ATi 9 £ 


SA6 


80+005010 


005176 








SA1 


80+005010 




SA2 


80+005011 


005177 








8X6 


^Xi*X2 




SA6 


80+005011 




NO 
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005200 
005201 



005202 

005203 

0D520A 

005205 
005206 



005207 
005210 



005211 



005212 



005213 



00521<* 



SA1 
SA2 

FX1 
NXO 
BX6 
LX6 

SA6 



SA1 

SA2 
1X1 
HXO 

BX6 
LX6 
SA6 



SA1 
SA2 

BX1 
MXO 
3X6 
NO 

SA6 



SA1 

BXO 
LXO 
BX1 
MXO 

BX6 
LX6 

SA6 



SA1 
LX1 
NO 

SA2 
LX2 
8X0 

LX2 
1X2 
1X6 
NO 



B0+005QQ** 
80 + 005005 

X2-X1 
01 

"■"XIaXO 
01 

B0 + QO5C1O 



60*005001 

B0 + 0Q50Q2 

X2-X1 

01 

"XIaXO 

01 

80+005010 



80+005010 
80+005011 

^X2*X1 
73 

"Xl-XO 



80+005011 



80+005001 

XI 
73 
Xl-XO 

01 

XOaXI 

01 

80+005011 



80+005001 
03 



80 + 005002 
01 

X2 



02 

X2+X0 

X1+X2 





SA1 


80+005001 


005216 








AX1 


03 




BXO 


xo-xo 




1X1 


xi+xo 




NO 




005217 








SA2 


80+005007 




AX2 


01 




1X6 


XI- X2 


005220 








SA6 


80+005002 




SA1 


80+005001 


005221 








3X0 


XI 




AXO 


Qk 




LXO 


0<* 




1X6 


XI- XO 


005222 








SA6 


80+005003 




SA1 


B0+OO5OO1 


005223 








SX6 


Xl+000144 


00522** 


SA6 


B 0+005007 








SA1 


B5+B0 




SB7 


Xl+80 




JP 


B7+000000 


005225 


17204000000000000000 


005226 


17216220771740 156064 



005215 



SA6 BQ+0Q5001 
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